Electronic Band gaps and transport properties inside graphene 
superlattices with one-dimensional periodic squared potentials 

Li-Gang Wang^'^ and Shi-Yao Zhu^'^'^ 
Centre of Optical Sciences and Department of Physics, 
The Chinese University of Hong Kong, Shatin, N. T., Hong Kong, China 
^Department of Physics, Zhejiang University, Hangzhou, 310027, China 
^Department of Physics, Hong Kong Baptist University, 
Kowloon Tong, Hong Kong, China 
(Dated: 2010.02.11) 

Abstract 

The electronic transport properties and band structures for the graphene-based one-dimensional 
(ID) superlattices with periodic squared potentials are investigated. It is found that a new Dirac 
point is formed, which is exactly located at the energy which corresponds to the zero (volume) 
averaged wavenumber inside the ID periodic potentials. The location of such a new Dirac point is 
robust against variations in the lattice constants, and it is only dependent on the ratio of potential 
widths. The zero-averaged wavenumber gap associated with the new Dirac point is insensitive to 
both the lattice constant and the structural disorder, and the defect mode in the zero-averaged 
wavenumber gap is weakly dependent on the insident angles of carriers. 

PACS numbers: 73.61.Wp, 73.20.At, 73.21.-b 
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I. INTRODUCTION 



Recently, the experimental realization of a stable single layer of carbon atoms densely 
packed in a honeycomb lattice has aroused considerable interest in study of their electronic 
properties Such kind of material is well known as graphene, and the low-energy 

charge carriers in pristine graphene are formally described by a massless Dirac equation with 
many unusual properties near the Dirac point where the valence band and conduction band 



■ dispersion, the chiral behaviour, ballistic 
6], frequency- dependent conductivity 



touch each other [l|-|5|], such as the linear eners 
conduction and unusual quantum Hall effect [3|, 
gate-tunable optical transitions j8|, and so on. 

Most recently, there have been a number of interesting theoretical invesitigations on the 
graphene supperlatices with periodic potential or barrier structures, which can be generated 
by diffierent methods such as electrostatic potentials js iSl and magnetic barriers 13-16|. 



Sometimes periodic arrays of corrugations 



17 



18| have also been proposed as graphene 



superlattices. It is well known that the supperlattices are very sucessful in controlling the 
electronic structures of many conventional semiconducting materials (e.g. see Ref. 19|). In 
graphene-based superlattices, researchers have found that a one-dimensional (ID) periodic- 
potential superlattice may result in the strong anisotropy for the low-energy charge carriers' 



group velocities that are reducedJ;o zero in one direction but are unchanged in another |10 |. 
Furthermore, Brey and Fertig [20] have shown that such behavior of the anisotropy is a 
precursor to the formation of further Dirac points in the electronic band structures and new 
zero energy states are controlled by the parameters of the periodic potentials. Meanwhile, 
Park et al. ^| pointed out that new massless Dirac fermions, which are absent in pristine 
graphene, could be generated when a slowly varying periodic potential is applied to graphene; 
and they further found the unusual properties of Landau levels and the quantum Hall effect 
near these new Dirac fermions, which are adjustable by the superlattic potential parameters 
22| . Finally it should also be mentioned that the electronic transmission and conductance 
through a graphene-based Kronig-Penney potential have been recently studied 23[ and the 
tunable band gap could be obtained in graphene with a noncentrosymmetric superlattice 
potential 24 1. 

Graphene superlattices are not only of theoretical interest, but also have been experi- 
mental realized. For example, superlattice patterns with periodicity as small as 5 nm have 
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been imprinted on graphene using the electron-beam induced deposition 25| . Epitaxially 
grown graphenes on metal (ruthenium or iridium) surfaces 26|-|32|| also show superlattice 
patterns with several nanometers (about 3 nm to 10 nm) lattice period. Fabrication of pe- 
riodically patterned gate electrodes is another possible way of making the graphene-based 
superlattices. 

Motived by these studies, in this paper, we will consider the robust properties of the 
electronic bandgap structures and transport properties for the graphene under the external 
periodic potentials by applying appropriate gate voltages. Following the previous work 
we evaluate the effects of the lattice constants, the angles of the incident charge carriers, the 
structural disorders and the defect potentials on the properties of electronic band structures 
and transmissions. It is found that a new Dirac point is exactly located at the energy 
with the zero (volume) averaged wavenumber inside the ID periodic potentials, and the 
location of such a new Dirac point is not dependent on lattice constants but dependent on 
the ratio of potential widths; and the position of the associated zero-averaged wavenumber 
gap near the new Dirac point is not only independent of lattice constants but is also weakly 
dependent on the incident angles. With the increasing of the lattice contants, the zero- 
averaged wavenumber gap will open and close oscillationaly but the center position of this 
gap does not depend on the lattice constants. Furthermore it is shown that the zero- 
averaged wavenumber gap is insensitive to the structural disorder, while the other opened 
gaps in ID periodic potentials are highly sensitive to the structural disorder. Finally we also 
find that the defect mode inside the zero-averaged wavenumber gap is weakly dependent on 
the insident angles while the defect mode in other gaps are highly dependent on the incident 
angles. 

The outline of this paper is the following. In Sec. II, with the help of the additional two- 
component basis, we introduce a new transfer matrix method, which is different from that in 
Ref. [ul, to calculate the reflection, transmission, and the evolution of the wave function; our 
transfer matrix method is very useful to deal with the periodic- or multi-squared potentials. 
In Sec. Ill, the various effects of the lattice constants, the incident angles of carriers, and the 
structural disorders on the electronic band structures are discussed in detail; furthermore 
the transport properties of the defect mode inside the zero-averaged wavenumber gap are 
also discussed. Finally, in Sec. IV, we summarize our results and draw our conclusions. 
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II. TRANSFER MATRIX METHOD FOR THE STRUCTURES OF PERIODIC 
POTENTIALS IN THE MONO-LAYER GRAPHENE 



The Hamiltonian of an electron moving inside a mono-layer graphene in the presence of 
the electrostatic potential V{x), which only depends on the coordinate x, is given by 



H = VF(T-p + V{x)I, (1) 

OX ' ay ' 



where p = {px^Py) = {~'^^'§^i~'^^'§~) is the momentum operator with two components, 



cr = ((7^,0"^), and ax,<^y are pauli matrices of the pseudospin, J is a 2 x 2 unit matrix, 
and vf ~ 10®m/s is the Fermi velocity. This Hamiltonian acts on a state expressed by a 
two-component pseudospinor \E' = {ip^^ip^Y , where ifjj^ and ip^ are the smooth enveloping 
functions for two triangular sublattices in graphene. Due to the translation invariance in 
the y direction, the wave functions 4'A,B{^:y) can be written as ^A,B{^iy) = i^A^Bi^)^^^^^- 
Therefore, from Eq. ([1]), we obtain 



kylpA = ikips, (2) 



dx 

+ kyi)B = ikipA, (3) 



dx 

where k = [E ~ V{xy\/hvF is the wavevector inside the potential V{x), E is the incident 
energy of a charge carrier, and k^ = E/Hvf corresponds to the incident wavevector. Obvi- 
ously, when E < V{x), the wavevector inside the barrier is opposite to the direction of the 
electron's velocity. This property leads to a Veselago lens in graphene p — n junctions, which 



33| 



has been predicted by Cheianov, Fal'Ko and Altshuler 

In what follows, we assume that the potential V{x) is comprised of periodic structures 
of squared potentials as shown in Fig. 1. Inside the j th potential, Vj{x) is a constant, 
therefore, from Eqs. ([2]) and (jS]), we can obtain 

+ (ij - kl)i., = 0, (4) 

^ + (*] - *J)*B = 0. (5) 

Here the subscript "j" denotes the quantities in the j th potential. The solutions of Eqs. 
(jl]) and ([5]) are the following forms 

^a{x) = ae"'^^ + 6e-'*", (6) 
^Ijb{x) = ce**-^ + (ie-^''^^ (7) 



where qj =sign{kj) ^/ kj — ky is the x component of the wavevector inside the j th potential 



Vj for kj > ky, otherwise qj = iyky — k^; and a (c) and b (d) are the amphtudes of the 
forward and backward propagating spinor components. Substituting Eqs. and ([7]) into 
Eqs. ([2]) and ([3]), we can find the relations 

a, (8) 



iqj + k, 



d = --^^b. (9) 

iqj ky 

Using Eqs. ([8]) and ([9]), we may obtain 

^^(x) = ae'*" + 6e-'*", (10) 

^^(a;) = a . e^-^^" - b . e"^'?^". (11) 

In order to derive the connection for the wave functions ipj^ ^{x) between any two positions 

Xj-i and Xj-i + Ax in the j th potential, we assume a basis $(a;) = (^^l^)), which are 
expressed as, 

cj)^(x) = ae'*" + fe-'"^", (12) 

(j)^{x) = ae'*" - fe-'"^". (13) 



Using the above basis, we can re-write Eqs. ( Iiup and fllip as the following form: 



where 



RjiE,ky)=\ . (15) 



i sin 9j cos 



Here =arcsin(A;y/A;j) is the angle between two components qj and ky in the j th potential. 
Inside the same potential, from the positions Xj-i to Xj^i + Ax, the wavefunction (^'^(!^^~|)) 
evolutes into another form {^^^^{~\'^^^-^) , which can be also expressed in terms of the above 
basis ^fx), 



V'a(xj_i + Ax)\ , .fM^j-i) 



where 



cos((7iAx) zsin((7,Ax) , , , 

T,(Ax,E,y=| ^ ^ |. (17) 

i sin(gj Ax + 9j) cos(gj Ax + 9j) 
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Therefore, the relation between ('^-^p and f^^)^-'' ^"J" .^0 can be finially written as: 



where the matrix Mj is given by 



Mj{Ax,E,ky) = Tj{Ax,E,ky)R-\E,k 



cos{qjAx—9j) ■sin{qjAx) 
cos 6j cos 6j 

■ sin{qjAx) cos{qjAx+9j) 



(19) 



It is easily to verify the equality: det[Mj] = 1. Here we would like to point out that in the 
case of = Vj , tl 
and it is given by 



case oi E = Vj^ the transfer materix (fT9|) should be recalculated with the similar method 



Mj{Ax,E,ky) = I ^ I . (20) 

exp{—kyAx) 

Meanwhile, in the j th potential < x < Xj), the wave functions ipy^ ^{x) can be also 

related with 4'a,b{^o) by 

where Axj = x—Xj^i, "ipABi^o) wave functions at the incident end of the whole structure, 
and the matrix Q is given by 

i-i 

QiAxj, E, ky) = Mj{Ax„ E, ky) JJ Mi{wi, E, ky). (22) 

i=l 

Here Wi is the width of the i th potential, and the matrix Q is related to the transformation 
of the charge particle's transport in the x direction. Thus we can know the wave functions 
^ab(^) position x inside each potential with the help of a transfer matrix. The 

initial two-component wave function (J^j^oj) can be determined by matching the boundary 
condition. As shown in Fig. 1, we assume that a free electron of energy E is incident from 
the region x < at any incident angle. In this region, the electronic wave function is a 
superposition of the incident and reflective wavepackets, so at the incident end (x = 0), we 
have the functions ^/'^(O) and ^/'b(O) follows: 

V'^(O) = Xlj,{E,ky)+,P,{E,ky) 

= {l + r)ij,{E,ky), (23) 



where ilj^{E,ky) is the incident wavepacket of the electron at x = 0. In order to obtain 
the function ipsi^) ^^e incident end, from Eqs. f|T2|) and f lT3|) . we can re- write the two- 
component basis in terms of the incident wavepacket, which are given by 

MO) = i^^{E, ky) + ME, ky) = (1 + r)ME, ky), (24) 
02(0) = ME, ky) - ME, ky) = (1 - r)ME, k.). (25) 

Therefore, using the relation of Eq. f|T4|) . we obtain 

^b(O) = i sin ^001 (0) + cos 0002(0), 

= {e'^'' -re~M4^i{E,ky), (26) 

where Oq is the incident angle of the electron inside the incident region (x < 0). In the above 
derivations, we have used the relation ip^{E,ky) = rip^^E.ky), where r is the reflection 
coefficient. Obiviously, we have 

^^(0)\^/ 1 + r 
In the similar way, at the exit end we have 



i^,{EM)- (27) 









V ^^(a^e) j 







ME,ky), (28) 

with the assumption of '0^(xe) = tilj^{E,ky), where t is the transmission coefficient of the 
electronic wave function through the whole structure, and dg is the exit angle at the exit 
end. Suppose that the matrix X connects the electronic wave function at the input end with 
Eq. (1271) and that at the exit end with ( l28l) . so that we can connect the input and output 
wave functions by the following equation 

xf^'-<»'V (29) 

\i'B(0) j 

with 

X= I ) = l[M^{w,,E,ky). (30) 

\ X21 X22 J j=i 

By substituting Eqs. fl27ll28p into Eq. fl29|) . we have the following relations 

t = {1 + r)xn + (e'^° - re-'^«)xi2, (31) 
te'^^ = + r)x2i + {e'^" - re-Mx22- (32) 
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Sovling the above two equations, we find the reflection and transmission coefficients given 
by 

""^ ' (x22e-^^o + xne^^O - Xue'('^-'^^ - X21 ' ^ ^ 

2 cos 6*0 .„,s 
' " (x22e-*^o + a:,,e^e^) - xi2e*(^-^") - X21 ' ^ ^ 

where we have used the property of det[X] = 1. With Eqs. (^3), ( 12^ and fl2^ . now we 
are able to calculate the two components of the electronic wave function at any position as 
follows 

^^(x) = i;,iE, ky) [(1 + r) Qn + (e^'° - re-^^«)gi2] , (35) 
^^(x) = V'.li?, ky) [(1 + ^) + (e^'" - re-*^«)Q22] , (36) 

where Qmn are elements of matrix Q. When we consider the translation of the electron in the 

y direction for obtaining the wave functions iPa,b{^i y)i the above functions i^A^Bi.^) have to 

i-i 

be producted by the factor, Y{y,ky) = exp (ikyAyj) Y[ exp(zfcj,Ayj), where Ay^ = WitanOi 

1=1 

and Ayj = Axj tanOj. Therefore we finally get 

^^(x, y) = ^,(E, ky)Y{y, ky) [{1 + r) Qn + (e^'" - re-'''')Qi2] , (37) 
^^{x, y) = ^.{E, ky)Y{y, ky) [(1 + r) Q^i + (e^'« - re~''^^)Q22] • (38) 

We emphasize that these equations are very useful to fully describe the evolutions of the 
two-component pseudospinor's wave function inside the potential or barrier structures of 
graphene when the incident electron's wavepacket ip^{E,ky) is given; and furthermore all 
these formulae are also suitable for multi-potential structures in graphene, not only for 
periodic potential structures. In the following discussions, we will discuss the properties 
of the electronic band structure and transmission for the graphene-based periodic squared 
potential structures. 



III. RESULTS AND DISCUSSIONS 



In this section, we would like to discuss some unique properties of the band structures in 
the graphene-based periodic-potential systems by using the above transfer method. First, 
let us invetigate the electron's bandgap for an infinite periodic structure (AB)^ , where the 
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periodic number N tends to infinity. The magnitude and width of the potential A (B) are 
with the electrostatic potential Va(b) and width wa{b), as shown in Fig. 1. According to 
the Bloch's theorem, the electronic dispersion at any incident angle follows the relation 

2cos[/3^D] = TtIMaMb] 

^ [2cOs{eA-0B)-2] , ^ , X . ^ 

= 2 cos[gAWA + qBWB] ^ 7 7 sm(gAWyi) ^vhKubWb)- (39) 

cos V A COS (7^ 

When the incident energy of the electron satisfies Vb < -E < Va, we have < 0, < 0, 
6'^ > 0, and > for the propagating modes. Then if qA^A = —(1bU!b, the above equation 
(|39|) becomes 

cos[/3,D] = 1 + ^-^^^^l^\MqAWA)\'. (40) 
cos^ Oa 

This equation indicates that, when qawa = —QbWb '^tt and 9a 7^ 0, there is no real 
solution for {3^, i.e., exsiting a bandgap; Note this bandgap will be close at normal incident 
case {9a = 0) from Eq. fj40|) . Therefore, the location of the touching point of the bands 
is exactly given by qAWA = —Qbuib at 9a = 0, i. e., kAWA = —kBUis or {Va — E)wa = 
{E - Vb)wb. 

Figure 2 shows clearly that a band gap opens exactly at energy E = 25meV under 
the inclined incident angles (i.e., ky 7^ 0), where the condition qAdA = —QBdB 7^ is 
satisfied. At the case of normal incidence {9a = 9b = 0), the upper and lower bands linearly 
touch together and form a new double-cone Dirac point. The location of the new Dirac 
point is governed by the equality: {Va — E)wa = {E — Vb)wb- For the graphene-based 
periodic-barrier structure with Va ^ and Vb = 0, the new Dirac point is exactly located 
at E = Va/(1 + wb/wa)- From Fig. 2(a-c), one can also find that the location of the new 
Dirac point is independent of the lattice constants; and the position of the the opened gap 
associated with the new Dirac point is not only independent of the lattice constants but also 
is weakly dependent on the incident angles [also see the transmission curves in Figs. 3a and 
3b for the finite structures, for example, {ABY^]] while other bandgap structuers are not 
only dependent on the lattice contants but also strongly dependent on different angles (i. e., 
different ky). The properties of the opened gap associated with the new Dirac point are very 
similar to that in the one-dimensional photonic crystals containing the left-haned materials 



34j |. where the so-called zero (volumn) averaged index gap is independent of the lattice 



constant but only dependent on the ratio of the thicknesses of the right- and left-handed 



materials 



34j |. From Fig. 4, one can also find that the locations of both the new Dirac 



point and the opened gap are determined by the ratio value of wa/wb for the cases with 
the fixed heights Va of the potentials. From the above discussion, we find that the volume- 
averaged wavenumber at the energy of the new Dirac point is zero, therefore such a opened 
gap associated with the new Dirac point may be called as the zero-averaged wavenumber 
gap. 

From Figs. 2(a-c) and 5(a-c), one can also noted that the slope of the band edges near 
the new Dirac point gradually becomes smaller as the lattice constant increases under the 
fixed ratio wa/wb and the fixed potential height, and such phenomena have been pointed 
out in recent work by Brey and Fertig 2o|. Actually, from Eq. ( 139|) . we can see that as the 
values of QaWa and —qs^B gradually reach the values of mvr (m = 1,2,3, ■■ ■), the slope 
of the band edges near the new Dirac point becomes smaller [see Figs. 2(a-c) and 5(a-c)]; 
once the condition qaWa = —(IbWb = "^tt is satisfied, the zero-averaged wavenumber gap 
will be close and a pair of new zero-averaged wavenumber states emerges from ky = [see 
Fig. 5(c)]. Here we would like to emphasize that the properties of these novel zero-averaged 
wavenumber states are similar to that of the zero-energy states in the previous work |20 |. 
Figure 5(d) shows how the zero-averaged wavenumber gap does gradually close with the 
increasing of the lattice constant under the fixed transversal wavenumber ky = O.Olnm^^ 
and the fixed ratio wa/wb- From Fig. 5(d), one can find that the zero-averaged wavenumber 
gap is very speical and it is open and close oscillationly with the increasing of the lattice 
constant, and the center position of the zero- averaged wavenumber gap is also independent of 
the lattice constant. However, the other opened gaps are largely shifted with the increasing 
of the lattice constant. 

Now we turn to consider the transmission of an electron passing through a finite graphene- 
based periodic-potential structure, e. g., {ABY'^, with the width deviation under different 
incident angles. Figure 6 shows the effect of the structural disorder on electronic transimitiv- 
ities. Figures 6a, 6b and 6c correspond to the transmissions through the periodic-potential 
structures with the width deviation (random uniform deviate) of ±2.5nm, ±3.75nm, and 
±5.0nm, respectively. From Figs. (6a)-(6c), it is clear that the higher opened gap is de- 
stroyed by strong disorder, but the zero-averaged wavenumber gap survives. The robustness 
of the zero-averaged wavenumber gap comes from the fact that the zero-averaged wavenum- 
ber solution remains invariant under disorder that is symmetric (+ and — equally probable). 
It should be emphasized again that the position of the zero-averaged wavenumber gap near 
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the new Dirac point is insensitive to both the incident angles and the disorder, [see Fig. 3b 
and Fig. 6]. 

Finally, we turn our attention to the effect of a defect potential on the property of 
the electron's transport inside the zero-averaged wavenumber gap. Here we consider the 
transmission of an electron passing through a graphene-based periodic potential structure 
with a defect potential, e. g., {AB)^^D{BA)^^, where the symbol "D" denotes the defect 
barrier. In Fig. 7(a), it shows the electronic transmitivity through the structure under 
different incident angles. One can see that there is a defect mode respectively occuring inside 
the zero-averaged wavenumber gap and the higher bandgap, and the location of the defect 
mode inside the zero-averaged wavenumber gap is very weakly dependent on the incident 
angle but the defect mode in the higher bandgap is strongly sensitive to the incident angle. 
It is clear shown in Fig. 7(b) that the energy of the defect mode inside the zero-averaged 
wavenumber gap is almost independent of the angles while the location of the defect mode 
in the higher bandgap has a large shift with the increasing of the incident angle. 



IV. CONCLUSIONS 



In summary, we studied the electronic band structures and transmission of carriers in the 
graphene-based ID periodic squared-potential superlattices. For the ID periodic squared- 
potential structure, it is found that a new Dirac point does exactly occurs at the energy that 
corresponds to the zero (volume) averaged wavenumber inside the system, and the location 
of such a new Dirac point is independent of lattice constants but dependent on the ratio of 
potential widths. It is also shown that the location of the zero-averaged wavenumber gap 
associated with the new Dirac point is not only independent of lattice constants but is also 
weakly dependent on the incident angles. As the lattice constant increases for the structures 
with the fixed ratio of potential widths, this zero-averaged wavenumber gap is open and close 
oscillationaly around the energy of the new Dirac point. We further showed the robustness 
of the zero- averaged wavenumber gap against the structural disorder, which has a sensitive 
effect on other opened gaps of the system. Finally we saw that the defect mode inside 
the zero-averaged wavenumber gap is weakly dependent on the insident angles while the 
defect mode in other gaps are highly dependent on the incident angles. Our analytical and 
numerical results on the properties of the new Dirac point, the novel band gap structure 
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and defect mode are hopefully of use to pertinent experiments. 
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Figures Captions: 

Fig. 1. (a) Schematic representation of the finite periodic squared potential structure 
in X — y plane. Dark regions denote the electrodes to apply the periodic potentials on the 
graphene. (b) The profiles of the periodic potentials applied on the monolayer graphene. 

Fig. 2. (Color online). Electronic band structures for (a) wa — wb — 20nm, (b) 
wa — Wb — 30nm and (c) wa — wb — 40nm, with Va — 50meV and Vb = in all cases. 
The dashed lines denote the "light cones" of the incident electrons, and the dot line denotes 
the location of the new Dirac points. 

Fig. 3. (Color onhne). Transmitivities of the finite periodic-potential structure {AB)"^^ 
under (a) different lattice contants with a fixed ratio wa/wb = 1 and an incident angle 
6q — 10° and (b) different incident angles with the fixed lattice parameters Wa = Wb = SOnm. 

Fig. 4. (Color online). Electronic band structures for (a) wa/wb — 1, (b) wa/wb — 3/2, 
and (c) wa/wb — 2, with Va — 50meV, Vb = and wb — 20nm in all cases. The dashed 
fines denote the locations of the new Dirac points. 

Fig. 5. (Color online). Electronic band structures for (a) wa = wb = 60nm, (b) wa = 
Wb = SOnm and (c) wa = wb = lOOnm; and (d) dependence of the bandgap structure on the 
lattice constant w with a fixed transversal wavenumber ky = O.Olnm"^ and wa = wb = w. 
Other parameters are the same as in Fig. 2. 

Fig. 6. (Color onhne). The effect of the structural disorder on electronic transimitivity 
T — under different incident angles: the solid lines for the incident angle 9o — 1°, the 
dashed fines for 9q = 5°,the dashed-dot lines for Oq = 10°, and the short-dashed fines for 
^0 = 15°, with Va = 50meV and Vb = 0, and wa = Wb = (20 + i?)nm, where i? is a random 
number for the case (a) between +2.5 nm and -2.5 nm, for the case (b) between -1-3.75 nm 
and -3.75 nm, and for the case (c) between -|-5 nm and -5 nm. 

Fig. 7. (Color onhne). (a) The transmitivity as a function of the incident electronic 
energy in a periodic-potential structure with a defect potential D, {ABY^D{BAy°. (b) 
Dependence of the defect modes on the incident angles. The parameters of the defect 
potential are wd = 70nm and Vd = 50meV, and other parameters are Va = 50meV, 
Vb = 0,and wa = wb = 20nm. 
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FIG. 2: (Color online) Electronic band structures for (a) = wb = 20nm, (b) wa = wb = 30nm 
and (c) WA = wb = 40nm, with Va = 50meV and = in all cases. The dashed lines denote 
the "light cones" of the incident electrons, and the dot line denotes the location of the new Dirac 
points. 
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FIG. 1: (Color online) (a) Schematic representation of the finite periodic squared potential structure 
m.x — y plane. Dark regions denote the electrodes to apply the periodic potentials on the graphene. 
(b) The profiles of the periodic potentials applied on the monolayer graphene. 
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FIG. 3: (Color online) Transmitivities of the finite periodic-potential structure (AB)"^^ under (a) 
different lattice contants with a fixed ratio wa/wb = 1 and an incident angle = 10° and (b) 
different incident angles with the fixed lattice parameters wa = wb = 30nm. 




FIG. 4: (Color online) Electronic band structures for (a) wa/wb = 1, (b) wa/wb = 3/2, and (c) 
wa/wb = 2, with Va = 50meV, Vb = and wb = 20nm in all cases. The dashed lines denote the 
locations of the new Dirac points. 
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FIG. 5: (Color online) Electronic band structures for (a) wa = wb = 60nm, (b) wa = wb = 80nm 
and (c) WA = wb = lOOnm; and (d) dependence of the bandgap structure on the lattice constant 
w with a fixed transversal wavenumber ky = O.Olnm^^ and wa = wb = w. Other parameters are 
the same as in Fig. 2. 
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FIG. 6: (Color online) The effect of the structural disorder on electronic transimitivity T = |tp 
under different incident angles: the solid lines for the incident angle Oq = 1°, the dashed lines 
for 6*0 = 5°, the dashed-dot lines for 6q = 10°, and the short-dashed lines for Oq = 15°, with 
Va = 50meV and Vb = 0, and wa = wb = (20 + i?)nm, where i? is a random number for the case 
(a) between +2.5 nm and -2.5 nm, for the case (b) between +3.75 nm and -3.75 nm, and for the 
case (c) between +5 nm and -5 nm. 
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FIG. 7: (Color online) (a) The transmitivity as a function of the incident electronic energy in a 
periodic-potential structure with a defect potential D, (AB)^^ D{B A)^^ . (b) Dependence of the 
defect modes on the incident angles. The parameters of the defect potential are wd = 70nm and 
Vd = 50meV, and other parameters are Va = 50meV, Vb = 0,and wa = wb = 20nm. 
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